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ABSTRACT 

We present detailed maps of the microlensing optical depth and event density over an 
area of 195 deg 2 towards the Galactic bulge. The maps are computed from synthetic 
stellar catalogues generated from the Besangon Galaxy Model, which comprises four 
stellar populations and a three-dimensional extinction map calibrated against the Two- 
Micron All-Sky Survey. The optical depth maps have a resolution of 15 arcminutes, 
corresponding to the angular resolution of the extinction map. We compute optical 
depth and event density maps for all resolved sources above I — 19, for unresolved 
(difference image) sources magnified above this limit, and for bright standard candle 
sources in the bulge. We show that the resulting optical depth contours are dominated 
by extinction effects, exhibiting fine structure in stark contrast to previous theoretical 
optical depth maps. Optical depth comparisons between Galactic models and opti- 
cal microlensing survey measurements cannot safely ignore extinction or assume it 
to be smooth. We show how the event distribution for hypothetical J and K-b&nd 
microlensing surveys, using existing ground-based facilities such as VISTA, UKIRT or 
CFHT, would be much less affected by extinction, especially in the K band. The near 
infrared provides a substantial sensitivity increase over current /-band surveys and a 
more faithful tracer of the underlying stellar distribution, something which upcoming 
variability surveys such as VVV will be able to exploit. Synthetic population models 
offer a promising way forward to fully exploit large microlensing datasets for Galactic 
structure studies. 
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^ 1 1 INTRODUCTION 



Several microlensing survey teams have been monitoring 
millions of stars over a l arge region of th e Galactic bulge 
for more than a decade (ISumi et al. 2003; Pop owski et all 
l2005l : ISumi et al.ll2006l : iHamadache et alJl2006h . These sur- 
veys have detected thousands of events and the combined 
detection rate of the OGLE-III and MOA-II surveys is cur- 
rently around 800 unique events per year. The microlensing 
optical depth, that is the instantaneous number of ongo- 
ing microlensing events per source star, is a key measurable 
for these surveys and provides an important constraint on 
the bulge surface mass density. Its dependence upon direc- 
tion provides, in principle, a powerful probe of the three- 
dimensional geometry of the bulge stellar mass distribution. 

To date several measurements have been made of the 
optical depth alo ng different directions. The MOA survey 
jSumi et al.ll2003h find an optical depth from difference im- 
age analysis (DIA) of 28 events of r di a = 2. 59^4 x 10~ 6 
for bulge and disc sources around Baade's Window (I ~ 
l°,b ~ — 3?9). Most recent measurements of the optical 



depth have been based upon a restricted subset of red 
clump giant (RCG) stars, which act approximately as stan- 
dard candles and are known to be confined to the bulge. 
From a sample of 42 RCG events the MACHO collabora- 
tion (|Popowski et al.ll2005l) measures a standard candle op- 
tical depth of T ac = 2.1 7+nlI x 10~ 6 at (I = 1?5, 6 = -2? 7). 
The OGLE-II survey (|Sumi et all 120061 ) reports a value of 
t sc = 2.55loie x 10 ~ 6 determined from 32 RCG events in 
the direction (I — l°2,b = —2? 7). The largest single sample 
analysed to date invo lves 120 RCG events ob tained by the 
EROS-2 experiment (|Hamadache et al"1l2006l h from which 
it measured the optical depth over a large region, includ- 
ing r sc = (2.42 ± 0.47) x 10 -6 within (-3?5 < b < -1?4). 
EROS-2 also reports an optical depth gradient with latitude 
|6] of dT sc /d\b\ ~ (0.78 ± 0.27) x 10~ 6 deg -1 at b = -2?7, in 
line with similar gradient determinations by MACHO and 
OGLE-II. 

Despite initial disagreement, optical depth measure- 
ments from RCG sources are now largely in accord 
with recent theoretical models (|Evans fc Belokurovl [20021 : 
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lHan fc Gouldll2003l : IWood fc Madl2005l ) and support the ex- 
istence of a bar-like bulge which is oriented 10 — 20° from 
the Sun — Galactic Centre line. However RCG events repre- 
sent less than 10% of all bulge microlensing events. A major 
limitation to the study of the remaining events comes from 
the current simplicity of theoretical models. A fully devel- 
oped model must account for the line of sight distribution 
of the sources in both the disc and bulge components. This 
distribution is a function not only of the underlying density 
models but also of the experimental sensitivity, the source 
luminosity function and the distribution of dust. The latter 
factor is particularly problematic for studies of the spatial 
distribution of events since extinction varies strongly with 
sky position. Until recently realistic three dimensional ex- 
tinction maps were unavailable. 

In this paper we present the first microlensing 
optical depth maps to incorporate a realistic three- 
dimensional extinctio n map. The map was developed by 
iMarshall et~"ai1 d2006h as part of the Besa ncon Galaxy Model 
IjRobin. Revle. Derriere fc Picaudl 120031 ) . In Section [5] we 
briefly overview the Besangon Model including how the 
three-dimensional dust map was generated. In Section [3] we 
describe how we set up the optical depth calculations. In 
Section [4] we present maps of optical depth and of event 
density for current /-band microlensing surveys, as well as 
for a hypothetical near-infrared survey which, as we show, 
would be much less affected by dust. We end with a discus- 
sion in Section [5] 



2 THE BESANQON GALAXY MODEL 

The Besangon Galaxy Model is a simulation tool aimed at 
testing galaxy evolution scenarios by comparing stellar dis- 
tributions predicted by these scenarios with observational 
constraints, such as photometric star counts and kinemat- 
ics. A comp lete description of the model ing r edien ts can 
be found in iRobin. Revle. Derriere fc Picaudf (|2003h . The 
model assumes that the Sun is at a Galactocentric distance 
Ro = 8.5 kpc. We summarise here the model's principal 
features which are relevant to the present analysis. In par- 
ticular, we do not discuss the details of the kinematics and 
metallicity of each population as the microlensing optical 
depth does not depend on them. 



2.1 Stellar populations 

The model assumes that stars are created from gas following 
a star formation history and an initial mass function; stel- 
lar evolution follows evolutionary tracks. To reproduce the 
overall galaxy formation and evolution we distinguish four 
populations of different ages and star formation history. 

• The spheroid popu lation is described in 



IRobin. Revle fc Crezel (|2000) as a single burst popula- 
tion of age 14 Gyr. Its density in the central regions is small 
compared to other populations and it is therefore only a 
marginal contributor to the optical depth. 

• The thick disc is formed from stars born about 11-12 
Gyr ago in a short period of time as implied by recent metal- 
licity determinations for this population. We also assume a 



single burst for simplicity. For the thick disc, star forma- 
tion occurred from the gas already settled in the disc. Kine- 
matics, deduced from observational constrain ts l|Qiha et al.l 
ll996l : IOiha. Bienavme. Mohan fc Robinlll999l '). imply that it 
has u ndergone a mergin g event shortly after the disc forma- 
tion ([Robin et al.l [l996). increasing the disc thickness and 
giving a higher velocity dispersion and scale height. The ini- 
tial mass function (IMF), density distributio n and local nor- 
malis ation are constrained from star counts (Rev le fc Robinl 
l200ll V 

• The bulge population is present in the centre of the 
Galaxy and extends to about 3-4 kpc. Its age is around 10±2 
Gyr. The population is modeled from iBruzual fc Charlotl 
(2003) population synthesis, as a single burst and a Salpeter 
IMF. This pop ulation reflects a t riaxia l bar distribution, as 
determined by IPicaud fc Robml (|2004j ) from near-infrared 
star counts. The angle of the major axis of the bulge with 
the sun-Galactic Centre axis is found to be 11?3. The scale 
lengths are 1.590 kpc, 0.424 kpc and 0.424 kpc. The total 
bulge mass is 2.05 x 10 10 M Q . 

• For the thin disc a standard evolutionary model is used 
to compute the disc population, comprising an IMF, a star 
formation rate (SFR) and a set of evo l utionar y tracks, as de- 
scribed in lHavwood. Robin, fc Crezel |l997bl ) and references 
therein. The disc population is assumed to evolve during 10 
Gyr, with multiple periods of star formation creating seven 
distinct age populations. A set of IMF slopes and SFRs are 
tentatively assumed and then tested against star counts. 
The tuning of disc parameters against relevant observa- 
tional d ata has been described in lHavwood. Robin, fc Crezel 
(|l997al lbh. A revised IMF has been adjusted to agree with 
the most recent Hi pparcos result s : the age-velocity disper- 
sion relation is fro mlGomez et al.l (Il997l). t he local luminos- 
ity function from Ijahreffi fc Wielenl (|l997h . giving an IMF 
slope a = 1.5 in the low mass r ange (0.5-0.08 Mq), in good 
agreement with iKroupal l|200lf ). The scale height has been 
computed self-consistently using the potential via the Boltz- 
mann equation. T he local dynamical mass is taken from 
ICreze et all l|l998h . 

Stellar remnants are also included in the model, though 
their role in the optical depth analysis is only indirect in 
contributing to the overall mass of lensing stars and in lim- 
iting the overall mass normalisation for visible stars which 
act as microlens sources; the microlensing optical depth is 
insensitive to the stellar or remnant mass function of the 
lenses (see Section |3j). We assume all white dwarfs to be of 
type DA and use evolutionary tracks and a tmosphere mod- 
els from_ BCTgeron^Teggettl^^uii (|200lh . complemented 
bv lChabrierl (|l999h for the very cool end, applicable to the 
halo. No brown dwarfs are included in the present model. 
All simulated stars are single stars. 

The evolutionary model fixes the distribution of stars 
in the space of intrinsic parameters: effective temperature, 
gravity, absolute magnitude, mass and age. These parame- 
ters are converted into colours in various systems through 
stellar atmosphere models. 

The detectability of stars computed in the simulations is 
subject to interstellar extinction and to observational errors. 
The latter are computed assuming Gaussian distributions 
with a dispersion as a function of the star magnitude. The 
extinction distribution is obtained as described below. 
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2.2 Extinction distribution 

iMarshall et al.1 (|2006h have built a 3D model of the distri- 
bution of the dust in the inner Galaxy (|Z| < 100°, \b\ < 10°) 
from an analysis of the 2MASS survey l)Skrutskie et al.l 
2006) using the Besangon Galaxy Model. By comparing the 
observed reddened stars to the unredenned simluated stars 
of the Besangon model, they were able to calculate the ex- 
tinction as a function of distance for any given line of sight in 
the Galaxy. The extinction along a line of sight is that which 
minimises the chi-squared difference between the J — K s 
colour distributions of the observed and simulated stars. The 
final map has an angular resolution of 15 arcmins and a dis- 
tance resolution of 100-500 pc. Better angular resolutions 
are possible but at the cost of reducing the resolution in 
distance. 

Having established the Galactic model we now turn our 
attention to its microlensing properties. 



3 OPTICAL DEPTH EVALUATION 

The optical depth is one of the primary measurables which 
can be obtained from statistical analyses of microlensing 
events. It depends only on the distribution and mass nor- 
malisation of the source and lens stellar populations and 
is therefore potentially a very powerful probe of Galactic 
structure. The o ptical depth to a source at distance S is 
l|Paczvriskilll986h 

T = [*-§ md L, (1) 

where L is the observer-lens distance, M is the lens mass, p 
is the lens mass density and Re = [4GML(S - L)/c 2 S] 1/2 
is the Einstein radius of the lens, with G and c having their 
usual meaning. 

An evaluation of r must involve an average over the 
observed source population. Such an average must take ac- 
count of, among other factors, the survey sensitivity. Cur- 
rent microlensing surveys such as OGLE-III and MOA-II 
are observing in the I band and are capable of 4% photom- 
etry at approximately J = 19 (c.f. Figure 5 of lUdalski et al.l 
|2002| ). We therefore adopt this limit for our /-band calcu- 
lations. We shall also consider a hypothetical near-infrared 
microlensing survey observing in the J or K band. We as- 
sume it to be capable of 4% photometry at J = 18 and 
K = 17 for comparison. Such a survey could be conducted 
with facilities such as UKIRT, VISTA or CFHT and, as we 
shall see, would have the advantage of being much less af- 
fected by dust than the /-band surveys. One survey which 
will be capable of detecting bulge microlensing events in the 
near-infrared is the Vista Variables in the Via Lactea (VVV) 
survey, which is scheduled to commence soon after the com- 
pletion of VISTA in 2009. 

We start by computing optical depths according to 
equation ([T]) over a three dimensional grid of source posi- 
tions (l,b,S), where I and b are Galactic coordinates. For 
the lenses we employ the same underlying mass density dis- 
tributions as are used in the Besangon Model simulations. 
We then use the Besangon population synthesis program to 
generate synthetic catalogues of millions of source stars. The 
catalogues cover the region \l\ < 9?75, |6| < 5°, giving a total 



area of 195 deg 2 . Each catalogue probes 3239 lines of sight 
to give a map resolution of 15' x 15', which is the limit of 
the resolution of the dust map. Finally, we interpolate the 
optical depth grid at the position of each synthetic source 
to arrive at a source- averaged optical depth, r, for each line 
of sight. 

We consider three different optical depth quantities. 
The first is r rcs , which is the optical depth averaged over all 
"resolved" sources, i.e. those which at baseline are brighter 
than the assumed survey limit of / = 19, J — 18 or 
K = 17. The second optical depth quantity is Tdi a which 
represents the average over all "difference image analysis" 
(DIA) sources, i.e. those which at peak magnification are 
brighter than the assumed survey limit. DIA sources can be 
significantly fainter than the survey limit at baseline and so 
our simulated source catalogues go four magnitudes fainter 
than our assumed survey limits. We weight the optical depth 
contribution of each DIA source by min(l,« 2 ), where u% is 
the largest impact parameter, in units of _Re, which allows 
the event to be detected above the survey limit. This is nec- 
essary because, at any instant, only a fraction u 2 of unre- 
solved sources are expected to be magnified sufficiently to 
be detectable as events. Many recent optical depth measure- 
ments have been based upon a subset of events involving 
bright red clump giant (RCG) sources, which are approxi- 
mate standard candles. We therefore compute a third opti- 
cal depth quantity, r sc , where r is averaged over a subset 
of bright "standard candle" sources. Our standard candle 
source selection has 

Io < min[/o, RC + 1, 9(V - /)o + JS.bc - 5.5], (2) 

where subscript denotes de-reddened magnitude and the 
red clump magnitude /o,rc — 14.66 — 0.033Z, with I the 
Galactic longitude expressed in degrees. This prescription 
follows the "exte nded red clump" selection algorithm of 
ISumi et al.1 (|2006l ). where we have adopted a weighted linear 
fit to the zero-point correct ed red clump mag nitude values, 
Johcj listed in Table 2 of ISumi et alJ l)2006h for Galactic 
latitudes b < 0. 



4 OPTICAL DEPTH AND EVENT DENSITY 
MAPS 

We present the results from Section[3]m two forms. Firstly as 
straightforward /-band maps of optical depth and secondly 
we present /-band maps of event density, which we define as 
the product of optical depth and source number density. The 
latter provides an instantaneous snapshot of the distribution 
of microlensing events on the sky. In the last part of this 
section we also present near-infrared microlensing maps. 

4.1 Maps of Optical Depth 

In Figure [l] we present /-band optical depth maps for re- 
solved sources (r rcs ), DIA sources (rdi a ) and bulge "standard 
candle" sources (r sc ). The optical depth is indicated both by 
the greyscale level and by the contours, which are set at 1, 
2 and 4 x 10" 6 . 

The most obvious distinguishing feature of these maps 
from previously published predictions is the irregularity of 
the contour line shapes. This is a direct consequence of the 
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l (deg) 

Figure 1. /-band optical depth greyscale maps, (a) The optical depth, r rcs , to all sources brighter then / = 19 at baseline, (b) The 
optical depth, T^ia, to all sources brighter than / = 19 at peak magnification, (c) The optical depth, t sc , to "standard candle" sources 
as defined by equation (0. Solid, dashed and dot-dashed contours indicate optical depths of 4, 2 and 1 X 10 — 6 , respectively. The low 
optical depth towards the Galactic plane is due to the high extinction in this region. The greyscale range is given by the bar to the right 
of each plot and is normalised to a maximum value of 7 X 10 -6 for all panels. 



effect of dust. At latitudes |6| < 1° the very high column den- 
sity of dust results in a negligible optical depth because only 
relatively nearby sources are detectable and these have an 
intrinsically small lensing optical depth from equation ([T]). 
The region of high extinction centred at I ~ —2°, b ~ —2° is 



also clearly evident from the microlensing optical depth dis- 
tribution. Also noticeable is the steep optical depth gradient 
towards the plane at \b\ ~ 1°. Away from the Galactic plane 
the optical depth con tours are smoother and comparable to 
previous models (e.g. lWood fc Maoll2005l ). In particular the 
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Figure 2. /-band grayscale maps of instantaneous event density, that is the product of optical depth and source number density (rn»). 
The units for rn* are events/deg 2 . Panels (a), (b) and (c) correspond to the cases plotted in Figure^ For comparison the optical depth 
contours in Figure [T] are also overplotted here in white. The top panel also shows 57 OGLE-III field locations (black squares) along with 
1069 events with Einstein radius crossing times 10 < is/days < 100 detected by the OGLE Early Warning System between 2005 and 
2007 (black dots). Note that the grayscale normalisation is different in each of the plots. 



predicted "standard candle" optical depth gradient south 
of the plane and away from regions of high dust density is 
dT BC /d\b\ ~ 0.7 x 10 -6 deg -1 , in ver y good agreement with 
measurements from RCG sources dPopowski et all 120051 ; 
ISumi et al . 2006; Hamadach e et al 1 l2006t l. 



In the absence of the dust the central hotspot in the 
optical depth and event density maps would be located at 
negative Galactic longitude due to th e higher optica l dept h 
towards the far side of the bar (e.g. IWood fc Mad [2005). 
However, this is much more difficult to see in the presence of 
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the dust distribution. When the dust is included the /-band 
optical depth peaks in a localised hotspot centred around 

I ~ 0?5, b ~ —2° with a value of 7x 10 -6 for standard candle 
sources in the bulge, 5.7 x 10~ 6 for resolved sources and 
6 x 10~ 6 for D IA sources. Evidence for such a hotspot comes 
from EROS-2 (|Hamadache et al.ll2006l) and from the O GLE- 

II Early Warning System fEWS. lUdalski et alj|2003h . The 
EROS-2 bright star event sample clearly shows the greatest 
contribution from the measured o ptical depth coming from 
events in this area (c.f. Figure 14 in lHamadache et al. 2006). 
In Section T4. 21 we discuss the distribution of events from the 
OGLE-III EWS. 

One interesting feature which emerges from Figure [T] is 
that the DIA optical depth Tdi a appears less prone to the 
effects of extinction, with a noticeably narrower avoidance 
zone due to dust near the Galactic plane. Also interesting are 
differences in the optical depth between the three measures 
for specific locations. For example, at Baade's Window (I — 
l°,b — — 3?9), both Tdia and r sc yield optical depths r ~ 
2.5 x 10" 6 whilst r ros ~ 1.5 x 10" 6 . 

Comparison of the op tical depth maps with recent opti- 
cal depth measurements llPopowski et all 120051 ; ISumi et al.l 
120061 ; lHamadache et all I2006T ) indicates a tendency of the 
model to overestimate the optical depth of bulge standard 
candle sources towards the plane at (at b > —3°). One 
factor which may be contributing to this is the relatively 
large Galactocentric distance scale employed in the model 
(Ro = 8.5 kpc). Equally, this tendency may also indicate 
that the bulge model and resulting extinction map may not 
be optimal. Another contributing factor may be that the 
extinction map effectively saturates in higher extinction re- 
gions (Av > 20) where the 2MASS isT-band sensitivity is 
relatively low. Consequently, the extinction may be under- 
estimated in regions where Ay > 15 (Schultheis et al, in 
preparation). Either way, a new bulge model for the Be- 
sangon simulation is in preparation and microlensing data 
should provide a useful consistency check for it. 

Another important factor when comparing predicted 
and observed optical depths is the photometric cuts used 
to define the source sample from which the optical depth is 
computed. The resolved source optical depth r rcs produces 
a significantly lower optical depth than either of the other 
two measures, highlighting the need for theoretical optical 
depth calculations to consider the photometric characteris- 
tics of individual survey selections. Whilst survey detection 
efficiency calculations may account for factors such as colour 
selections, they nonetheless assume some underlying lumi- 
nosity function for the microlensing sources, which may not 
necessarily be the same as assumed in theoretical calcula- 
tions. 

4.2 Maps of Event Density 

Figure [2] shows equivalent /-band maps for the instanta- 
neous event number density, which we define as the prod- 
uct of source-averaged optical depth and the source number 
density. This quantity is indicative of the spatial variation 
in the event rate provided there are not significant spatial 
gradients in the average event duration. The grayscale map 
in the plots in Figure [2] trace the event density and, to aid 
comparison, we have overplotted the optical depth contours 
from Figure [T] In the upper panel of Figure [5] we plot the 



positions of 57 OGLE-III fields together with 1069 events 
detected within these fields between 2005 and 2007 by the 
OGLE-III EWS3. We have excluded events with Einstein 
radius crossing times £e < 10 days in order to suppress the 
bias induced by the fact that some of the OGLE-III fields 
are monitored several times per night whilst others are mon- 
itored nightly. Additionally we exclude long-duration events 
with £e > 100 days so as to minimise the risk of tracing 
kinematic structures that would not be represented in our 
maps which are based only on optical depth. 

For all the /-band maps there is a predicted hot-spot 
in the event density at around I ~ l°,fe ~ —2°. The maps 
clearly indicate that one expects many more events south of 
the plane than north of it, supporting the current strategy of 
the survey teams in focussing predominately on the southern 
side of the plane. 

Interestingly the distribution of OGLE-III EWS alerted 
events also shows hotspots in the regions predicted by the 
model, and generally the model correlates quite well with the 
alert event distribution. Of course the detection efficiency 
of the EWS events is completely uncalibrated so one must 
be wary of drawing too much from such a comparison, but 
on the face of it the model appears capable of reproduc- 
ing clumpy features seen in the observed event distribution. 
One exception may be at I ~ 1?5, b ~ — 1?5 where the model 
predicts a strong /-band microlensing signature for resolved 
sources but which is not covered by the OGLE survey, which 
has been optimised over 16 seasons in order to concentrate 
on high event yield locations. Nonetheless, with sufficiently 
sophisticated models there is clearly the exciting prospect 
of using the large available samples of microlensing events 
to constrain fine structure in the underlying Galactic struc- 
ture. However, to do this the detection efficiency of the event 
samples must be known. 

4.3 Near-infrared microlensing 

We finish by considering the potential i mpact o f mi- 
crolensing surveys in the near- infrared |Gouldl 1 19951 : 
lEvans fc BelokurovS 2002 ) . Specifically, we consider the event 
density in the J and K -bands for resolved sources, assum- 
ing a survey capable of 4% photometry at J = 18 and 
K = 17. Such a survey could be performed using the 
wide-field infrared cameras available on telescopes such as 
UKIRT, VISTA or CFHT. Indeed the Vista Variables in 
the Via Lactea (VVV) survey is scheduled to commence on 
VISTA in 2009 as an ESO Public Survey. VVV will survey 
some 520 deg 2 of the disk and bulge, primarily in the K 
band, over five seasons, with at least one season of almost 
nightly monitoring from which a substantial number of mi- 
crolensing events should be detected. 

The top, middle and bottom panels in Figure [3] show, 
respectively, the /-band map (also shown in the top panel 
of Figure [5}, the J-band map and the Zf-band map of event 
density for resolved sources (i.e. using r ros ). The grayscale 
levels of the maps are matched for direct comparison and, as 
in Figure [5J the contours indicate the optical depth. The ad- 
vantage of observing in the near-infrared is very clear. Firstly 
the effects of extinction are greatly reduced, especially in K, 

1 http : //www . astrouw . edu . pl/~ogle/ogle3/ews/ews . html 
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so that the event density is much higher and, as attested by 
the smoothness of the AT-band optical depth contours, more 
faithfully traces the shape of the underlying mass distribu- 
tion. These two factors alone mean that a /f-band survey 
like VVV would have powerful advantages over /-band sur- 
veys for Galactic structure studies. However the other main 
advantage of observing in J or K is that the northern hemi- 
sphere also becomes an attractive target for microlensing 
studies. In particular the region around |Z| < 2°, 6 ~ 1?5 
exhibits a strong microlensing hotspot in the K band, and 
a reasonably strong hotspot in J. The significance of this is 
that this region includes the far side of the bar at negative I. 
The ability to survey both the near and far sides of the bar 
would allow the separate microlensing contributions of the 
disc and bulge to be separated, as the far bar should have a 
significantly higher disc optical depth than the near bar. 



5 DISCUSSION 

We have used a synthetic population synthesis model of the 
Galaxy to construct maps of the microlensing optical depth 
and event density. The Besangon Model employs a three- 
dimensional extinction map and allows optical depth maps 
to be computed for individual passbands, and for the maps 
to adopt sophisticated cuts in source colour and magnitude. 

We have computed /-band optical depths maps for three 
distinct event samples: stars which at baseline are resolved 
at / < 19 (T r es), comparable to current microlensing survey 
limits, stars which reach / < 19 at peak magnification (rdia), 
and sources which are selected using colour cuts similar to 
those used to define standard candle samples such as red 
clump giants (r sc ). The structure in the resulting optical 
depth maps is dominated by the effects of extinction to- 
wards the Galactic pane, and demonstrates a much stronger 
microlensing signal South of the Galactic plane, which is 
where most of the survey effort is being concentrated. We 
find that rdia tends to be less affected by the dust. More 
generally, the comparison between the observed and mod- 
elled optical depth can be quite sensitive to the photometric 
characteristics of the event selection. 

We also compute maps of event density, that is the prod- 
uct of optical depth and source density. We find good corre- 
lation between these maps and the location of 1069 events 
detected by the OGLE-III Early Warning System (EWS). 
A direct quantitative comparison cannot be made since the 
efficiency of the EWS is unknown. However the model suc- 
cessfully predicts the region of highest event yield as seen by 
both the OGLE-III and EROS-2 surveys, suggesting that 
the large scale features of the model are correct. We also 
find good agreement between the predicted optical depth 
gradient with latitude and that determined by microlensing 
surveys. However, the model presented here tends to over- 
predict the bulge optical depth towards the plane in the case 
of standard candle source samples. This may possibly be due 
to a combination of a relatively large assumed Galactocen- 
tric distance and an underestimation of extinction in high 
extinction regions. An improved Galactic model is currently 
being developed. 

We also investigate the potential of near-infrared mi- 
crolensing surveys by computing event density maps in the J 
and K bands. The K band in particular is extremely promis- 
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Figure 3. A comparison of /, J and /C-band maps of instan- 
taneous event density. The grayscale in each map is normalised 
to a peak density of 15 events/deg 2 . Panel (a) is a reproduc- 
tion of panel (a) in Figure [2] whilst panels (b) and (c) shows the 
event density in the J and K bands, respectively, for all sources 
brighter than J = 18 and K = 17. Such a near-infrared microlens- 
ing survey could be undertaken with UKIRT, VISTA or CFHT. 
The advantage of lesser extinction, particularly in the K band, 
is clearly evident. The contours indicate optical depth (not event 
density), with levels set the same as for Figure [l] 

ing in circumventing the worst effects of interstellar dust. K- 
band optical depth contours more faithfully trace the under- 
lying structure of the disc and bulge and also allow both the 
near- and far-bar to be probed for microlensing. Comparison 
of near- and far-bar microlensing samples would allow the 
separate microlensing signatures of disc and bulge lenses to 
be decoupled, providing additional constraints on Galactic 
models. Upcoming near-infrared variability surveys, such as 
VVV, have the potential to provide excellent microlensing 
samples for Galactic structure studies. 

The number of discovered microlensing events is ap- 
proaching 5000, and yet so far only a few hundred have been 
analysed for Galactic structure studies. Whilst theoretical 
models which ignore extinction (or assume it to be smooth) 
can be compared reasonably successfully with samples of 
~ 100 events, they likely will be inadequate for dealing with 
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samples of several thousand events, where fine structure may 
start to become evident. An accurate calibration of the stel- 
lar mass budget in the inner Galaxy is also a vital prerequi- 
site to evaluating the allowable mass budget of dark matter 
in the inner Galaxy. The analysis of microlensing datasets 
can therefore provide important indirect constraints on the 
dark matter halo profile. We believe that synthetic modelling 
provides an excellent way forward for such work. By fully ex- 
ploiting the large and expanding catalogues of microlensing 
samples, synthetic modelling provides the potential for mi- 
crolensing to be used as a high-precision probe of Galactic 
structure. 
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